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Abstract 

Binary systems comprising at least one neutron star contain strong grav- 
itational field regions and thereby provide a testing ground for strong-field 
gravity. Two types of data can be used to test the law of gravity in com- 
pact binaries: binary pulsar observations, or forthcoming gravitational- wave 
observations of inspiralling binaries. We compare the probing power of these 
two types of observations within a generic two-parameter family of tensor- 
scalar gravitational theories. Our analysis generalizes previous work (by us) 
on binary-pulsar tests by using a sample of realistic equations of state for 
nuclear matter (instead of a polytrope), and goes beyond a previous study 
(by CM. Will) of gravitational-wave tests by considering more general tensor- 
scalar theories than the one-parameter Jordan-Fierz-Brans-Dicke one. Finite- 
size effects in tensor-scalar gravity are also discussed. 

PACS numbers: 04.80.Cc, 04.30.-w, 97.60.Gb 
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I. INTRODUCTION 



The detection of gravitational waves by kilometric-size laser-interferometer systems such 
as LIGO in the US and VIRGO in Europe will initiate a new era in astronomy. One of the 
most promising sources of gravitational waves is the inspiralling compact binary, a binary 
system made of neutron stars or black holes whose orbit decays under gravitational radiation 
reaction. The observation of these systems will provide important astrophysical information, 
e.g. masses of neutron stars, and direct distance measurements up to hundreds of Mpc U . It 
is also said that detecting gravitational waves from inspiralling binaries should provide rich 
tests of the law of relativistic gravity in situations comprising strong-field regions (like near 
a neutron star or a black hole). However, present binary-pulsar experiments already provide 
us with deep and accurate tests of strong-field gravity [§H|J. ^ * s therefore interesting to 
compare and contrast the probing power of present (and foreseeable) pulsar tests with that 
of future gravity-wave observations. 

A convenient quantitative way of doing this comparison is to work within a multi- 
parameter family of physically motivated (and physically consistent) theories of gravity 
which differ from Einstein's theory in their radiative and strong-field predictions. The most 
natural framework of this type is the general class of tensor-scalar theories in which gravity 
is mediated both by a tensor field g^ u ("Einstein metric") and by a scalar field ip. These 
theories contain one arbitrary "coupling function" A((p) defining the conformal factor relat- 
ing the pure spin-2 Einstein metric g* to the "physical metric" g^ u = A 2 (ip)g^ u measured 
by laboratory clocks and rods. The usual Jordan-Fierz-Brans-Dicke theory [|5],§],[7|] is the one 
parameter class of theories defined by a coupling function A((p) = exp(a v ? )- [The coupling 
strength a is related to the often used parameter u by = (2uj + 3) -1 .] Will || has stud- 
ied the quantitative constraints on the coupling parameter oq of Jordan-Fierz-Brans-Dicke 
theories that could be brought by gravitational- wave observations. His result is that in most 
cases the bounds coming from gravity-wave observations will be comparable to presently 
known bounds coming from solar-system experiments (namely, «q < 10~ 3 ). This result of 
Ref. || seems to suggest that gravitational-wave-based strong-field tests of gravity do not go 
really beyond the solar-system weak-field tests of gravity. We wish, however, to emphasize 
that this seemingly pessimistic conclusion is mainly due to having restricted one's attention 
to the special, one-parameter Jordan-Fierz-Brans-Dicke theory. Indeed, in this theory the 
strength of the coupling of the scalar field <p to matter is given by the constant quantity cto 
independently of the state of condensation of the gravitational source. As a consequence, 
the predictions of the theory differ from those of Einstein's theory by a fraction of order 
«q in all situations: weak- field ones or strong- field ones, alike. By contrast, it has been 
emphasized in Refs. PJTOll that the more generic tensor-scalar theories in which the strength 



of the coupling of ip to matter, namely 

^sf 1 ' (L1) 

depended on the value of ip, allowed for the existence of genuine strong-field effects, by which 
the presence of a highly condensed source, such as a neutron star, could generate order- 
unity deviations from general relativity even if the weak-field limit of a(ip) is arbitrarily 
small. [More precisely, these non-perturbative strong-field effects take place when @((p) 
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da((p)/dip is negative.] This led us to consider, instead of the Jordan-Fierz-Brans-Dicke 
model A((p) = exp(a ip), the class of theories defined by 

,%) = expQ/V) . (1.2) 

This class of theories contains two arbitrary (dimensionless) parameters: /3q appearing 
in Eq. (|1.2|), and ip , the asymptotic value of ip at spatial infinity. [By contrast, in Jordan- 
Fierz-Brans-Dicke theory, ipo has no observable effects.] Equivalently, the two parameters in 
these theories can be defined as being the strength of the linear coupling of if to matter in 
the weak-field limit (ip ~ <p ), 

dlnA(<p Q ) 

a = a{tpo) = — = (J Lp , (1.3) 

and the non-linear coupling parameter 

da(ip ) d 2 In A(<p ) 

PO — — o — o • \ LA ) 



A convenient feature of the two-parameter family of theories (|1.2|) is that they have 
just the amount of generality needed both to parametrize the most general boost-invariant 
weak-field ("post-Newtonian") deviations from Einstein's theory, and to encompass nonper- 
turbative strong-field effects. Indeed, on the one hand, the two theory-parameters (a ,/3 ) 
determine the two well-known Eddington-Nordtvedt-Will parameters, 

7 = 7Eddin g ton - 1 = -2 ajj/(l + al) , (1.5) 

P = Aldington - 1 = g A) "o/l 1 + a l) 2 > C 1 - 6 ) 

which measure the most general, boost-invariant, deviations from general relativity at the 
first post-Newtonian level. (See |TT| , |12"1| for the generalization of the Eddington parameters 
to the second post-Newtonian level.) On the other hand, when (3$ < —4 non-perturbative 
strong-field effects develop in the theories defined by Eq. ( |1.2|) . 

All existing gravitational experiments can be interpreted as constraints on the two- 
dimensional space of theories defined by Eq. (|1.2j) . In other words, we can work within 
the common (aoj A)) plane, and consider each gravitational experiment (be it of weak-field 
or strong-field nature) as defining a certain exclusion plot within that plane. In some recent 
work [0, we have constructed such exclusion plots, as defined by considering both solar- 



system experiments and binary-pulsar experiments. The present work will generalize these 
exclusion plots in several respects: (i) we shall plot the regions of the (a , /3 ) plane probed by 
future gravitational-wave observations of neutron star-neutron star, and neutron star-black 
hole systems, (ii) we shall improve our previous study of the probing power of binary-pulsar 
experiments by considering a sample of realistic nuclear equations of states (instead of the 
polytropic one used by us before) and by using updated pulsar data, and (iii) we shall con- 
sider the individual constraints on a and (3q brought by the main solar-system experiments 
(instead of using published combined limits on j and 0). 
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This paper is organized as follows: In Sec. II we summarize our (numerical) approach 
to computing the various form factors that describe the coupling of the scalar field ip to 
a neutron star. In Sec. Ill we generalize Ref. || in discussing how gravitational wave 
observations can give quantitative tests of tensor-scalar gravity. In Sec. IV we combine 
and compare gravitational-wave tests with binary-pulsar tests and solar-system ones. Our 
conclusions are presented in Sec. V, while an Appendix discusses finite-size effects in tensor- 
scalar gravity. 



II. GRAVITATIONAL FORM FACTORS OF NEUTRON STARS DESCRIBED BY 

REALISTIC EQUATIONS OF STATE 

The orbital dynamics of a binary system depends, besides the Einstein masses of the two 
objects rriA, rriB, on the effective coupling constants ct^, defined as 

d In m a 

a A = ^~ -, 2.1 

as well as on their scalar-field derivatives 

f>A=%±. (2.2) 

The derivatives in Eqs. ( |2.1| ), ( |2.2| ) are taken for fixed values of the baryonic mass mA- In 
the limit of negligible self-gravity a a — » «o an d Pa — > Po, but it was shown in Refs. PJIQ 



that (when p is negative enough) the effective coupling constant of a neutron star can 
reach values of order unity even if the weak-field coupling constant «o is extremely small. 
By contrast, in Jordan-Fierz-Brans-Dicke theory a a tends to zero with ocq. Let us also note 
that for a black hole = because of the no-scalar-hair theorems. 



In Ref. JTTJ it was shown that the parameters tua, m B , a a, o>b, Pa, Pb suffice to deter- 
mine both the (v /c) 2 -accurate conservative orbital dynamics (including periastron precession 
and other relativistic deformations of Keplerian motion) and the radiation reaction effects 
(monopolar, dipolar and quadrupolar effects linked to the combined emission of scalar and 
tensor waves). On the other hand, the relativistic timing of binary pulsar systems involves, 
besides the above effective coupling constants, a new parameter describing the possible field 



dependence of the (Einstein-frame) inertia moment I a of the pulsar [|13|,|TIJ. This new pa- 
rameter, entering the timing formula, is of the form K% = kA ol b where 

k A = 5 • (2.3) 

For a fixed value of the baryonic mass m A , the quantities m A , a a, Pa and kA depend 
on: (i) the value (fo of the scalar field at infinity and the theory of gravity used (hence a 
dependence on a and Pq in our case), and (ii) the equation of state of nuclear matter used 
to describe the interior of the neutron star. 

In this work, we shall consider, besides the simple polytropic equation of state that we 
used before, 
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~~~ K n a rrih f n\ ~ _ / n\ , 
e = nm b + — , P = Kn m b — , 2.4 

r — 1 Vr2 / \n / 



m fc = 1.66 x 10" 24 g , n = 0.1 fm" J , T = 2.34 , K = 0.0195 , (2.5) 

a sample of three (more or less) realistic equations of state. They describe the crust by 
the same equation of state (due to Baym et al. |]14|| ) but differ in their descriptions of the 
core of the neutron star. One description (due to Pandharipande ||15|| ) gives a rather soft 
equation of state, a second one (due to Wiringa et al. |TTJ) is of intermediate stiffness, while 
the third one (due to Haensel et al. |L7|]) corresponds to a somewhat stiff equation of state. 
This sample will allow us to see the effect of the softness of the equation of state. 

We have discussed in detail in Ref. |T(J the way to compute the effective coupling con- 



stants rriA, &a, Pa, in tensor-scalar gravity. We shall not repeat this material here. Let 



us only mention a modification that we brought to the procedure discussed in Ref. [10] for 
more convenience when dealing with equations of state given in tabular form. Starting from 
such tables, we first interpolate them to define the (proper) physical energy density e and 
the (proper) physical baryon number density n as functions of the physical pressure p. [The 
tilde refers here to quantities expressed in "physical units", i.e., with respect to the metric 
g^y = A 2 g^ u .] Then, we transform the field equations by using the pressure p as a radial 
variable. More precisely, this means that we consider a set of 8 first-order differential equa- 
tions for the variation with p of the variables M, v, ip, if), p, M, uj and w. The symbols 



M, . . . ,tu have the same meaning as in Ref. | 10| . In particular, p is the Schwarzschild-like 



radial coordinate used in Eq. (3.1) of |TIJ. Note that the quantity a is constructed at the 
end of the numerical work from the asymptotic value of ip by using a = (3q<Po, Eq. ( |1.3|) . 



The right-hand side of the equation giving dp /dp is the inverse of the right-hand side of 
Eq. (3.6e) of [Kj, while the right-hand sides of the 7 other equations is simply obtained 
[from dQ/dp = (dQ/ dp) /(dp/ dp)] by dividing the right-hand sides of Eqs. (3.6a)-(3.6h) by 
the right-hand side of Eq. (3.6e). A bonus of this modified approach is that we can numer- 
ically integrate these differential equations on a known range: namely, we start integration 
at the center of the star with some given central value for p, say p c , and we stop at p = 
(surface of the star). The value of p at the end of integration gives the radius R of the star. 
We have discussed in Refs. [p|,p!0| how to obtain the quantities m^, m^, a a, and I a from the 
other results of the integration. 

Actually, having chosen some nuclear equation of state, the above numerical integration 
yields, for each value of /3 , and for each initial values p c (central pressure) and ip c (central 
value of the scalar field), an output which consists of (m^, tha, a a, I a, <Po), where y?o is 
the value of ip at spatial infinity. [We keep only the "positively scalar-polarized" configura- 
tions with <po<p c > 0. Indeed, we have shown in ]10|] that the solutions with (po<p c < are 



energetically disfavored.] Having constructed a sufficiently dense grid (in (p c , (p c ) space) of 
such numerically integrated models, we can then compute Pa, Eq- ( |2.2|) , and Ica, Eq. Q2.3|) , 
by numerically approximating the y? -derivatives (with fixed m^) by finite differences con- 
structed from four neighboring grid points. The expression (|2.1| ) for a a is used as a test of 
our numerical accuracy (we appropriately densify our grid so as to ensure that Eq. ( [2.1| ) is 
satisfied to better than 5%). Using finally some interpolation, we end up by generating a 
huge, multi-entry table which is conveniently organized as giving, for a grid of values of «o 
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and Po, numerical values of (mA, ^A, oca, Pa, I a, &a)- This is a numerical approximation to 
giving rriA, a a, Pa, I a and Ua as functions of ac , p and rf%A- We should also mention that 
we keep only stable models in our tables, i.e., we keep only configurations corresponding to 
masses increasing, when the central pressure increases, between (formally) zero and their 
first maximum (corresponding to the usual concept of the maximum mass of a neutron star, 
here quantitatively modified by scalar- field effects). 



III. TESTING TENSOR-SCALAR GRAVITY FROM GRAVITATIONAL- WAVE 
OBSERVATIONS OF INSPIRALLING BINARIES 

The gravitational wave signal from inspiralling binaries will be deeply buried in the 
broadband noise of interferometric detectors. To detect it, one will have to correlate the 
output of the detector with a "template wave form", i.e., an hopefully faithful copy of the 
actual time evolution of the signal. As hundreds to tens of thousands of cycles may be 
usefully registered in the output of the detector, this method of matched filtering will allow 
one to dig deeply into the broadband noise, under the condition that the used template wave 
form be an accurate representation of the signal over that many cycles. This means that 
gravitational wave observations will be very sensitive to the evolution of the frequency and 
the phase of the wave. 

In tensor-scalar gravity the evolution of the phase of the gravitational wave will be 
different from that in Einstein's pure tensor gravity because this phase evolution is driven 
by radiation reaction, which is modified because the system can now loose energy to scalar 
waves, in addition to loosing energy to tensor waves. More precisely, the energy lost to 
scalar or tensor radiation at infinity now reads, with sufficient approximation, 

£j ^monopole _|_ ^idipole _|_ ^iquadrupole _|_ ^nquadrupole (3 1) 

where the subscript denotes the field being radiated, and where the superscript refers to the 
multipolar structure of the emitted radiation. The explicit expressions for the various energy 
losses have been derived for general tensor- (multi-)scalar gravity theories in Ref. |l l| (see 
also Hl8|). As said above, they depend on the (Einstein) masses uia, m B and on the effective 



scalar coupling parameters of the two stars, ola, ot B) Pa, Pb- Denoting M = itla + m B , 
Xa = ttia/M, X b = m B /M = 1— X A , v = XaX b , and assuming that the orbit is 
circularized (eccentricity e = 0), we find from Eq. (6.52) of Ref. |TT| the following evolution 
law for the orbital frequency n = 2n/P or ut 



n v 



n 2 1 + olaoib 




(3.2) 



where 



v = [G*(l + a A a B ) Mn] 1/3 (3.3) 



(in which G* denotes the "bare" gravitational constant entering the action of the Einstein 
metric: fd 4 xg¥ 2 R(g*)/16ir G*), and 
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re = 1 + - {a A X B + a B X A f + \ (a A - a B )(a A X A + a B X B )(X A - X B ) 
6 6 

+ 7^ T~7 — {PB<y A X A - (3 A a B X B ) + d 2 (a A - a B ) 2 . (3.4) 

48 1 + a A a B 

The term (a A — a B ) 2 {v/c) 3 in Eq. (|3~"2|) is associated to the emission of scalar dipolar 
waves |T3| , p!9| , p!8| , pO| . The first term (equal to one) in re comes from the emission of tensor 
quadrupolar waves, the second term is due to the emission of scalar quadrupolar waves, the 
remaining terms come from 0((v/c) 5 ) contributions to the scalar dipolar waves [responsible 
for the 0((v/c) 3 ) term in Eq. (|3.2| )1. The numerical coefficient d 2 has not been computed, 
but this is not important in the present context as it multiplies the same term (a A — a B ) 2 
which will be more strongly constrained by the 0((v/c) 3 ) dipolar radiation term in Eq. ( |3.2| ). 
[Let us note in passing that our parameter re differs by a certain factor from the one used 
by Will ||, namely, re wm = (1 + a A a B ) 2 (1 + ag)~ 3 re here .] 

The expression ( |3~2l) is similar to the frequency evolution equations used by Will ||. Our 
results are, however, more general because they apply in a generic tensor-gravity theory, 
while his were restricted to the Jordan-Fierz-Brans-Dicke case. On the other hand, starting 
from Eq. ( |3.2| ) we can follow his analysis in discussing the constraints that can be derived 
from (future) gravitational wave observations. Let us define a tensor-scalar "chirp mass" (in 
time units) by 

M = ^(l + a A a B )M (- ) , (3.5) 

and the corresponding adimensionalized orbital frequency 

u = Mn. (3.6) 
In terms of these variables, Eq. (|3.2|) reads 

M U = y « 11/3 { 1 + b V 2 ' 5 U 2 ' 3 + 0(b U 1 ' 3 ) 

_/743 11 \ 2/5 2/ 3[1 + ( a 2^ 
V336 4 J 

+4tt u~ 3/5 u[1 + 0(a 2 )} + 0(m 4/3 )|, (3.7) 



where we defined 

b = -| re- 3 / 5 (1 + a A a B )- 2 "> (a A - a B ) 2 . (3.8) 
96 

The zeroth-order term oc u 11 ^ 3 in Eq. fl3.7| ) represent the frequency evolution due to the 
emission of tensor quadrupolar waves. The other terms represent corrections to this zeroth- 
order result. The terms containing a factor b or a symbolic factor a 2 represent the effects 
of the radiation of scalar waves. The main term due to the radiation of scalar waves is the 
b z/ 2 / 5 u~ 2 ^ 3 term. It is important to note that this term (coming from dipolar waves) contains 
the large factor u~ 2 l 3 ~ v~ 2 ^> 1. This term is the main source of the constraints on tensor- 
scalar gravity coming from gravity-wave observations because it has a different dependence 
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on u than the usual general relativistic terms. Indeed, the usual general relativistic frequency 
evolution is given by the terms without b or a 2 factors, i.e., symbolically, M.ii ~ w 11 / 3 [1 + 

u 2/3 + w 3/3 + M 4/3 + ...]_ 

Assuming that future gravity-wave observations from inspiralling binaries will be well 
matched by filters constructed from the standard general relativistic orbital phase evolution, 
we wish to quantify what constraints they could bring on the modifications to the frequency 
evolution equation ( |3.7| ) by the terms proportional to b or a 2 . This problem has been tackled 
by Will H using a matched-filter analysis. His analysis can be applied to our case, if we 
follow his (plausible) assumption that the final bounds will restrict scalar-field effects to be 
more or less uniformly small compared to general relativistic ones. This means that the 
only essential scalar-field term that we should consider in Eq. (|3.7| ) is the term b z/ 2 / 5 w~ 2 / 3 . 
All the scalar-field effects that simply modify usual einsteinian effects can be neglected. 
Furthermore, this means also that we can approximate k by one, and the chirp mass ( |3.5|) 
by its einsteinian limit c~ 3 M z/ 3 / 5 . Finally, we can translate Will's final bound on the 
Jordan-Fierz-Brans-Dicke coupling parameter a 2 , (his Eq. (|1.5|) ), in the following bound on 
b or (a a — o>b) 2 '- 

« 6 „ I (aA - aB f < 1.46 x 10- ( £ ) ( " ) , (3.9) 



Aitational time associated to a solar mass. 



where T = G* m c is the gravitational time associated to a solar mass. This bound can 
also be rewritten in terms of M = itla + itlb and v = rriAfnB/i^A + it^b) 2 '- 



( M \ 7/3 / 10 \ 
(<*-„,)•< 2.92x10- - (3.10) 



In Eqs. ( p.9| ) and (|3.10| ), S/N denotes the signal-to-noise ratio (after applying a matched 
filter) of a LIGO/VIRGO-type detector. [The analysis of Ref. f| assumes the noise spectral 
density of a LIGO "advanced detector".] 

We have focussed above on the indirect effect of scalar-wave emission on the phase 
evolution of tensor waves for the following reason. Relatively to the main tensor-quadrupolar 
contribution, this effect is of order (a^ — as) 2 (c/v) 2 ^> (a^ — as) 2 - Moreover, because of 
possible nonperturbative scalar-field effects associated to the strong self-gravity of neutron 
stars (ctA— &b) 2 can reach values of order unity, even if the weak-field scalar coupling «o 1. 
By contrast, the more direct effect of scalar- wave emission on the differential forces acting 
on any (local) detector of gravitational waves can be described by the transverse projection 
of the physical metric g^ v 

% = (A 2 (cp + 0) « KJ T + 2 a 81 , (3.11) 

where Sfj = 5ij —rii rij is the projection of 5^ transversally to the direction n % of propagation 
of the wave (see, e.g., |]TTJ), and where we choose units such that A((p$) = 1. The main direct 
effect of scalar-wave emission, associated, to scalar dipolar emission, is then (see Eq. (6.33) 
of 0) 

2 a di ? olc 8% = -~ a (a A - a B ) n t {z\ - 4) . (3.12) 
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Compared to the standard tensor-quadrupolar wave, this is of order ao (ola — ocb) (c/v) 2 . 
Even when nonperturbative strong- field effects generate a a ~ &b = 0(1), we see that the 
direct effects of scalar-wave emission are down by a factor «o <C 1 compared to the indirect 
ones. One therefore a priori expects that these direct effects will be less sensitive probes 
of strong-field deviations from general relativity than the indirect ones considered above. 
This qualitative-argument-based conclusion is confirmed by the results of Will (Sec. IV of 
||) who quantified the lower sensitivity of matched filtering to amplitudes, as compared to 
phases. 

Let us also note that Refs. ||ID] have emphasized that gravitational collapse and neutron- 



star binary coalescence (by contrast with inspiral) might exhibit observably significant direct 
and indirect scalar-wave effects due to the combined possibilities of strong-scalar-field effects 
and monopolar radiation. [However, the recent investigation of Ref. |2l] suggests that these 
effects will be too small to be of observational interest.! 



IV. COMBINING GRAVITY- WAVE, PULSAR AND SOLAR-SYSTEM TESTS 

Thanks to our consideration of a two-dimensional "mini-space" of alternative gravity 
theories, we can combine various experimental tests of relativistic gravity, and compare 
their probing power. Indeed, each experimental constraint separates the («o,/3o) plane into 
two regions: an allowed one and an excluded one (at some fixed confidence level). We can 
therefore compare the probing power of various experimental tests by drawing a combined 
exclusion plot in the (a , /3 ) plane. 

We shall consider three solar-system tests of weak-field gravity. The measurement of the 



Shapiro time delay by the Viking mission |22| as well a some VLBI measurements p3|j24 
yield the lcr bound |t~| < 2 x 10~ 3 , which (in view of Eq. (|1.5|) ) gives 



a 2 < KT 3 (la) . (4.1) 



The measurement of the perihelion advance of Mercury |25] gives 3 \2j — (3\ < 10 3 which, 
using Eqs. (|1.5|) and Ql.6|) , translates into 



2 1 /O 2 



< 7.5 x KT 4 (la) . (4.2) 



Finally, the Lunar Laser Ranging experiment [EH] yields —1.7 xlO 3 < 4 /3 — 7<3xl0 4 
which translates into 

-8.5 x KT 4 < al + (3 Q al < 1.5 x KT 4 (la). (4.3) 

These three constraints, taken together, exclude the region above the solid line labeled 
"1PN" in Figs. [!]-§]. Concerning binary-pulsar experiments we shall take into account three 
of them. 

The PSR 1913+16 experiment uses the values of three well-measured phenomenological 
timing observables ( "post-Keplerian parameters" 0]) ci; obs , 7° bs and P b obs . [We neglect the 
low-precision measurement of the post-Keplerian parameters r and s] Here u> ohs denotes 
the observed secular rate of advance of the periastron, 7° bs denotes the observed value of 
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the "Einstein" time-dilation parameter (entering = 7sinw ||]), and P b obs denotes the 
secular change of the orbital period. The values we shall take for these observed parameters 
are (27| 

u ohs = 4.226 621 (11)° yr" 1 , (4.4a) 
7 obs = 4.295 (2) x 1(T 3 s , (4.4b) 
pobs _ pgai = _ 2 .4ioi(85) x 10~ 12 , (4.4c) 

where figures in parentheses represent (realistic) la uncertainties in the last quoted digits. 
In Eq. ( |4.4cj ), P b gal represents the sum of various galactic effects which must be subtracted 
from P b ohs to be able to compare it to theoretical predictions ||28|| . [We will take into account 



below the small modifications brought to it by tensor-scalar gravity, Eq. (9.22) of We 
need also the values of the Keplerian parameters 

P b = 27 906.980 780 4(6) s , (4.5a) 

e = 0.617 130 8(4). (4.5b) 



As explained in detail in Ref. |10] one converts the three phenomenological (theory- 
independent) measurements (|4.4j) into constraints on alternative gravity theories by us- 
ing the predictions that, say, tensor-scalar gravity make for u, 7, Pf, as functions of 
the two (a priori unknown) masses m^, mg and of the parameters (here ocq, Po) 
defining each tensor-scalar theory. As already mentioned above, the dependence on 
the theory-parameters aco, Po goes through the dependence of the post-Keplerian pa- 
rameters u), 7 and Pb on the effective coupling constants a a, oc b , Pa, Pb, k>A- The 
explicit expressions of uj theoiy (mA, tub, oca, Pa, ^a, ■ ■ ■), l theorj ( m A, m B,otA, Pa,^a, ■ ■ ■) and 
P b themy (rriA, tub, oca, Pa,^a, ■ ■ ■) have been written down in Refs. [|TTJ and [fLOR . We shall not 



reproduce them here. Using these expressions, we then define a goodness-of-fit statistics by 
minimizing over the physically a priori undetermined masses: 

X 2 (pf s ; oca, A)) = min mAjms \ 2 (pf s ; ™a, m B , oc A , Pa, k A , ■ ■ ■) , (4.6) 

where, for the simultaneous measurement of any number of post-Keplerian parameters pi, 
one would define 

X 2 ( P f s ; m A , m B , a A , Pa, ■ ■ ■) = E (Pi^i™* ■ ■ " P^Y ■ ( 4 ' 7 ) 

i 

Here, we neglect the correlations between the measurements of the various post-Keplerian 
parameters. 

In principle, when considering one or more pulsar experiments simultaneously, one can 
associate a confidence level (in the bayesian sense) to any region 1Z in the theory-plane by 
integrating over TZ the normalized a posteriori probability density 

da dpo 7r post (a , Po) = da dp N 7r prior (a , Po) Wi(p° bs | a , Po) ^(pf 8 | a , Po) ■ ■ ■ (4.8) 

where the probability (density) Wi to observe, in the first pulsar experiment, the multi- 
plet p° bs (e.g. (u ohs , 7° bs , P b ohs )) of post-Keplerian parameters can be approximated by the 
Gaussian 
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WMpf s K,A))=exp 



-~ x 2 (pf s ;« ,/3 ; 



(4.9) 



with the definition (|4.6|) above. In Eq. ( |4.8| ), AMs a normalization constant (such that 
/ J da dj3 vr post (a , Po) = 1) and 7r prior (ao, Po) is some a priori probability density for the 
values of the theory parameters. A common (and convenient) way of choosing some confi- 
dence regions TZ is to define them from the various x 2 ■ One can consider either the over- 
lap (assumed to exist) of the various individual x 2 contour levels x 2 (p° bs ; a o, Po) — ki, 
X 2 (P2 bs j a o, Po) = ^2, etc., or the combined x 2 contour levels: Xtot( a o, Po) = k to t, with 
Xtot = X 2 (Pa bs j a o, Po)- To each choice of 71, we associate a bayesian probability ("proba- 
bility of causes" ; the causes being here the unknown theory parameters a , /3 ) by integrating 



(4.8) over 7Z. Here, we shall not try to quantify confidence levels in this way. We shall content 
ourselves by plotting individual x 2 — 1 contours. If several independent so-defined pulsar 
contours overlap, it is clear that the overlap region will define a rather high confidence level 
region for the total experiment combining these independent pulsar measurements. 

As discussed in Section II above, it is essential in this discussion to decide upon some 
equation of state for nuclear matter. Each choice of nuclear equation of state defines its own 
corresponding exclusion plot. For example, we expect soft equations of state (like Pand- 
haripande's) to yield stronger constraints on the theory parameters («o> AO- Indeed, they 
lead to more highly condensed neutron star models, and thereby to generically higher values 
of the effective couplings a a, Pa, kA for given (a , Po) [jTTH - We shall consider successively 



the following equations of state: (i) in Fig. [I] the polytrope (|2 . 4j) - (P75|) (as used in our pre 



vious work |H|), (ii) in Fig. ^| the (soft) Pandharipande equation of state, (iii) in Fig. |3| 



the (medium) Wiringa et al. equation of state, and (iv) in Fig. ^ the (stiff) Haensel et al. 
equation of state. 

The PSR 1534+12 experiment consists mainly of the high-precision measurements of 
three post-Keplerian parameters (which have to deal with strong-field effects without mixing 
of radiative ones 0): cj obs , 7° bs and s obs . Here s obs denotes a phenomenological parameter 
measuring the shape of the gravitational time delay P5,4]|. The values we shall take for these 



three observable parameters are |30 



c^ obs = 1.755 76(4)° yr~\ (4.10a) 

7 obs = 2.066(10) x 10" 3 s, (4.10b) 

s ohs = 0.982(7) . (4.10c) 

We shall also take into account the lower-precision measurement of the range-of-time-delay 
parameter r obs 

r obs = 6.7(1.3) x 10" 6 s, (4.11) 

and we shall need the Keplerian observables 

P b = 36 351.702 6587(35) s , (4.12a) 

e = 0.273 677 7(5), (4.12b) 

x = 3.729 463(3) s. (4.12c) 
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Reference |30| gives also the observed value of the radiative parameter P b obs . However, 
it underlines that the theoretical significance of this measurement is (at present) highly 
uncertain because of the lack of an independent, reliable measurement of the distance to 



PSR 1534+12. Indeed, as discussed in pSfl , several galactic effects have to be subtracted from 
P b obs before equating it to the theoretical prediction p fe theory . These corrections are relatively 
small and sufficiently well known in the case of PSR 1913+16, while they are relatively large 
and insufficiently known in the case of PSR 1534+12. Because of this insufficient knowledge 
of the theoretically relevant combination P 6 obs — -P b gal , we shall not take into account this 
observable in our exclusion plots below. [An alternative choice would be to take it into 
account but to enlarge the errors induced by the uncertainty on the distance. In that case, 
we found that its main effect is to forbid the horn-shaped region at the top-left of Fig. [l], 
which is anyway already ruled out by solar-system experiments as well as other binary-pulsar 
tests. 1 



Finally, we follow Ref. |TI| in taking also into account the data on PSR 0655+64. This 
binary pulsar is composed of a neutron star of mass rriA ~ 1.4% and a white dwarf 
companion of mass ~ 0.8 m Q , moving around each other on a nearly circular orbit in 
a period of about one day. This dissymmetric system is, potentially, a strong emitter of 
scalar dipolar waves. Indeed, we saw above that dipolar radiation losses are proportional 
to («a — «b) 2 (f/c) 3 . Here «b does not differ significantly from the weak-field coupling ao 
because the self-gravity of a white dwarf is very small, while a a can reach values of order 



unity. We refer to Ref. [10] for the discussion of the constraints obtainable from this system. 



The end result is that we get the following (conservative) la constraint 

[a A (m A ) - « ] 2 < 3 x 10~ 4 . (4.13) 

We use the corresponding Xo655+64 = 1 level in our exclusion plots. 

For all the equations of state that we consider, we see on the exclusion plots Fig. [I]-41] that 
the most stringent constraints coming from pulsar experiments are obtained by combining 
the 1913+16 exclusion region for /3 < +1 with that from 0655+64 for j3 > +1. The 
resulting theoretical bounds are somewhat less stringent (by a factor of a few) than solar- 
system bounds when /? ^ /3 C , while for fi < /3 C pulsar experiments essentially exclude an 
infinite domain of the (ao, /3q) plane which remained allowed by solar-system experiments. 
Here, j3 c denotes the (negative) critical value of /3q below which nonperturbative strong- 
field effects develop, thereby exhibiting the unique strong-field probing power of pulsar 
experiments. As we see on Figs. [l|-f|, the value of f3 c depends on the equation of state. In 
particular, a soft equation of state leads to highly condensed neutron star configurations 
and, thereby, develop nonperturbative effects earlier than stiff equations of state: in other 
words — /^ oft < — (3f" lS . This is visible on Fig. ^| (Pandharipande) where the pulsar bound is 
more stringent than the solar-system ones for (3 Q < —1, while for stiffer equation of state it 
becomes more stringent only when fio < —3. 

Finally, we have added, for comparison, the exclusion regions defined by the gravity- 
wave observation limit ( |3.10|) , assuming a signal-to- noise ratio S/N = 10. In absence of 



gravity-wave observations telling us about the precise masses of real inspiralling binaries, 
we have considered two fiducial cases: (i) a two-neutron-star system with Einstein masses 
rriA = 1.441 m Q , ins = 1.388 m (as measured in PSR 1913+16 when interpreted in general 
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relativity), and (ii) a neutron star-black hole system with itia = 1.4 m , = 10 m®. In 
case (i) neither the precise numerical values of the masses nor the fact that we fix Einstein 
masses instead of baryonic masses is crucial. What is crucial in our definition of the fiducial 
case (i) is that we assume a fractional mass difference Am/m m Am/m ~ 4%, as large as in 
PSR 1913+16. Indeed, the over important dipolar radiation is proportional to [^(m^) — 
«b(^b)] 2 oc (Am/m) 2 as Am/m — > 0. In case (ii), the no-scalar-hair theorems guarantee 
that gib = for a black hole (see, e.g., [0]), so that neutron-star-black- hole systems are 



always good a priori probes of possible scalar dipolar radiation. 

Because of the complexity of the numerical calculation of the strong-field form factors of 
neutron stars (see Sec. II above), we could compute more precisely the exclusion regions for 
the polytropic equation of state, Eqs. Q2.4D - fl2.5D . Therefore, it is in Fig. [l] that one sees best 
the shape of the regions possibly excluded^ by future gravitational wave data. The fiducial 
case (i) (a la 1913+16) excludes an ellipsoidal bubble which touches the /3q axis around 
/3 Q « —5, while the fiducial case (ii) (neutron star-black hole) excludes the region above the 
nearly straight line a — 0.03 (3 ~ 0.15 represented in Fig. |l|. The corresponding excluded 
regions for realistic equations of state can be recognized on Figs. |2|-[| as deformations of the 
just described regions for the polytropic case. Note that the bubble excluded by 1913+16-like 
systems (i) is smaller when the EOS is softer, and that it even disappears for Pandharipande's 
equation of state (Fig. Q). In that case, the detection of such a system by a gravitational- wave 
interferometer would not constrain at all the space of theories. 

The value S/N = 10 chosen for Figs. |1}J4] corresponds to the conventional event rate of 3 
binary-neutron-star coalescences per year in a radius of 200 Mpc, and to a probably smaller 



event rate for neutron-star-black-hole coalescences. However, as pointed out by Will pi 



the event rate is only slightly relevant to our discussion since a single system can suffice to 
constrain the space of gravity theories. It is therefore interesting to consider also the lucky 
discovery of an exceptionally near system, with a signal-to-noise ratio as large as S/N = 100. 
The corresponding exclusion plots are displayed in Fig. |5], for the same polytropic equation 
of state as in Fig. |l|. The bubble excluded by the 1913+16-like system is much larger, but 
still not competitive with present binary-pulsar tests. Similarly, the neutron-star-black-hole 
system is more constraining than in Fig. [1], but the slope of the corresponding line is only 
reduced by a factor ~ V4~0, although the signal-to-noise ratio is 10 times larger. This is 
due to the fact that the dominant dipolar radiation is proportional to the square of a^. 
Therefore, the lucky discovery of a nearby neutron-star-black-hole system would be slightly 
more constraining than present binary-pulsar tests in the region — 3 ^ fio < 0, but not 
better than present solar-system experiments. Moreover, one must keep in mind that solar- 
system experiments will improve in the mean time. In particular, NASA's Gravity Probe B 
mission (due for launch in 2000) is expected to improve the probing of ao down to the level 
a ~ VlO" 5 w 3 x lO" 3 . 

Since we are mentioning the possibility of lucky discoveries, like those of PSRs 1913+16 
and 1534+12, let us also quote the constraints which could be achieved if a binary pulsar 



1 We assume here that general relativity is the (nearly) exact description of gravity chosen by 
nature, and we discuss constraints on deviations away from general relativity. 
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with a black-hole companion were discovered. In that case, the main dipolar radiation 
would be proportional to a 2 A (since «b = for a black hole), instead of the small factor 
(a a — &b) 2 oc (Am/m) 2 appearing for binary-neutron-star systems. The mass of the black 
hole is not a crucial parameter for this discussion (we take mg = 10 m Q ). Assuming an 
orbital period Pb and a measurement accuracy for Pf, similar to those of PSR 1913+16, one 
finds that the constraints on cyq would be tightened by a factor! ~ 80, i.e., that the bold 
lines of Figs. |T]-^] would cross the vertical axis around cto ~ 0.05/80 ~ 6 x 10 -4 . In terms 
of the Eddington parameter 7Edd, Eq. (|1.5| ), this corresponds to a level |7Edd _ 1| ~ 10~ 6 , 
which is about a thousand times tighter than present solar-system limits, and ten times 
better than the probing level expected from Gravity Probe B. This underlines that, in the 
future, binary-pulsar tests may become competitive with, or even supersede, solar-system 
experiments even in the region /3q > —4 of the theory plane. 



V. CONCLUSIONS 

The main conclusion of the comparison carried out in Figs. []]-§] is that, in all cases, future 
LIGO/ VIRGO observations of inspiralling compact binaries turn out not to be competitive 
with present binary-pulsar tests in their discriminating probing of the strong-field, and 
radiative, aspects of relativistic gravity. This conclusion may seem paradoxical. It should 
not be interpreted negatively against LIGO/VIRGO observations which, as shown in Figs. |l|- 
|5], will independently probe strong-field gravity and will exclude regions of parameter space 
allowed by solar-system experiments. Rather, it is simply a reminder that binary-pulsar 
experiments are superb tools for probing strong-field and radiative aspects of gravity. It 
is also somewhat a good news for gravitational wave data analysis (which promises to be 
already a very challenging task even if one a priori assumes the validity of general relativity; 
see, e.g., [j32|). Indeed, our results Figs. [L|-|5| indicate that our present experimental knowledge 
of the law of relativistic gravity is sufficient to justify using general relativity as the standard 
theory of gravitational radiation. 

Note that this conclusion explicitly refers to the quantitative probing of plausible! devi- 
ations from general relativity. At the qualitative level, and also at the non-discriminating 
quantitative level, LIGO/VIRGO observations will bring invaluable advances in our experi- 
mental knowledge of relativistic gravity. First, they will provide the first direct observation 
of gravitational waves far in the wave zone (while binary pulsar experiments prove the reality 



2 This improvement factor of 80 on gq comes mainly from the ratio |ons — obh|/|oa — o#| ~ 
|cns — cbh|/|ca — cb\, where the "compactness" parameters [11 for neutron stars and black holes 



are respectively ca ~ O.21m^/m0 and cbh = 1- (Here = 1.441 m , = 1.388 m , ttt-ns ~ 
1.4 m©, and the value of meH does not matter.) 

3 Within the presently accepted framework for low-energy fundamental physics, namely field the- 
ory, the only alternative (long-range) gravity theories which, (i) do not violate the basic tenets 
of field theory, and (ii) are not already necessarily extremely constrained by existing equivalence- 
principle tests, are the tensor-scalar gravity theories. 
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of the propagation with finite velocity of the gravitational interaction in the near zone of 
a binary system). Second, they will (hopefully) lead to superb additional confirmations of 
general relativity through the observation of the wave forms emitted during the inspiral and 
coalescence of neutron stars or black holes; see, e.g., [fll,j33l|. 

Independently of the comparison between LIGO/VIRGO tests and binary-pulsar tests, 
the present paper has provided the first systematic study of the influence of the nuclear 
equation of state on the theoretical probing power of binary-pulsar tests. In particular, Fig. |2| 
shows that if the equation of state were as soft as predicted by the simple Pandharipande 
model, binary-pulsar tests quantitatively supersede solar-system ones in all the region (3 < 
— 1 of parameter space. Even if we consider the less constraining stiff equations of states, 
the present work confirms the limit 

Po > -4.5 (5.1) 

found (modulo 10%) in [[K]]. We recall that this limit can be interpreted as a limit on the 
ratio of the two weak-field post-Einstein parameters 

£^«-i/3b<l.l. (5.2) 
7Edd - 1 4 

Finally, we pointed out that the discovery of a binary pulsar with a black-hole companion 
has the potential of providing a superb new probe of relativistic gravity. The discriminating 
power of this probe might supersede all its present and foreseeable competitors in measuring 
ao down to the level o& < 10~ 6 . 
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APPENDIX A: FINITE-SIZE EFFECTS IN TENSOR-SCALAR GRAVITY 

In the text we have assumed that the leading modification, in tensor-scalar gravity, of 
the orbital motion of binary systems comes from the change in radiation reaction forces. 
In this appendix, we briefly discuss the modifications of the orbital motion caused by the 
finite extension of the bodies. Contrary to the pure spin 2 theory where such effects are 
very small (because they are suppressed in spherical bodies), the presence of a scalar field 
opens the possibility of couplings to the spherical inertia moment / d 3 xp(x)x 2 . In this 
appendix we use the general diagrammatic approach of Ref. [|T^j to confirm (and hopefully 



better understand) the results of Nordtvedt Q.^.^g] based on some explicit calculations 
valid only for weakly- self-gravitating extended bodies. The final outcome is that finite-size 
effects can be neglected in a matched-filter analysis of inspiralling compact binaries. 
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Following the indications given in ||12|| , we can formally, but very generally, take extension 
effects into account by considering that the effective action for each compact body is a 
functional of the fields ip and g^ u which can be expanded in terms of the values along some 
central worldline of their spacetime derivatives (derivative expansion). Namely, we take the 
following action for the extended body (labeled E) 

^tended = _ Q J m f g*^ (-^(4) dz% dz^ , (Al) 

with 

™>f [<P, 9%] = m E ((p(z)) + I E (<p) R* + Je(<p) Ue u e R% + k e(<p) d * V 

+L E (<p) u% < V; d„<p + M E (<p) u*£ < d„ <p d vV + N E (<p) g*^ d^d u¥ >. (A2) 

This is the most general action, expanded up to two derivatives of the fields (p and/or 
for a body which is spherically symmetric and static when unperturbed. [Terms which are 
first order in derivatives are excluded by spherical symmetry or, for H E (<p) u* E <9 M (p, by time- 
reversal symmetry] This form, being generic, is valid for strongly self-gravitating bodies. 
In this case, the (^-dependent quantities I E , J E , K E , L E , M E , N E , define some "scalar form 
factors" of body E which go beyond the basic effective coupling a E ((p) = dlnm E ((p)/d<p 
associated with the point-like effective action S E omt = —cj m E ((p) d s* E . [For non-spherical 
and/or non-static bodies many other new scalar form factors could, in principle, appear.] 

Most of the a priori independent-looking terms in Eq. (|A2j) can be easily shown either not 
to contribute at the (observationally most relevant) first post-Keplerian level [i.e., 0(v 2 /c 2 ) 
beyond the Keplerian orbital motion], or to be equivalent to other terms, modulo some 
redefinition of the dynamical variables. Let us first recall that any correction term, in a 
Lagrangian, which is proportional to the zeroth-order equations of motion can be redefined 
away by shifting some of the dynamical variables. Technically, S [if>] + el (if)) 8S / 'Sip = 
So[ip'} + 0(e 2 ) with ip' = ip-\-eI(ip). In our case, the dynamical variables can be either g^ u , ip 
or z E . In all cases, the local redefinitions ip — ► ip' have no effects on the observables at infinity 
(such as the periastron advance or the evolution of the orbital phase of an inspiralling binary). 
Using such local redefinitions of the dynamical observables (and formally neglecting singular 
self-action terms, i.e., 5-function contributions to rn E [(p]), one easily checks the following 
simplifications: The Je(<^) contribution can be transformed (using the g^ u field equations 
R*^ v = 28^ p> d u ip + (^-functions and a local redefinition of g*^ u within body E) into the M E (ip) 
contribution. The M E (tp) term contributes only at the second post-Keplerian level 0(v 4 /c A ). 
The K E (tp) term can be eliminated by redefining ip within body E. The L E (ip) term can 
be reabsorbed in the N E (ip) one by redefining the position z E , and by neglecting second 
post-Keplerian (2PK) effects. Indeed, integrating by parts — / cds* E L E (ip) u* fJ, u* u V* d u ip 
one finds an Me(<^)-like term (contributing only at 2PK order) plus a term proportional to 
u *ti y * u * ^ j s e q Ua i to — a E (5^ + u* M u*) <9 M <p upon using the equations of motion. This 

shows that after the redefinition 

z E ^4 = z E + ^- V*<p , (A3) 
m E 

the L E (ip) contribution is, modulo 2PK, equivalent to changing N E (ip) — > N' E (<p) = N E + 
a E L E . Similarly, the I E term can be absorbed in the N E one by locally redefining g*. [In 
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fact, if one fixes the harmonic gauge before making this redefinition of g* , one also needs 
a shift of the positions of the bodies A ^ E to absorb I E into Ne, namely, z\ — > z'\ = 
4 + 4a/ £ (4-4)/(r| A c 2 )+0(l/c 4 ).] 

Finally, we end up (modulo 2PK) with a generic, simplified effective action containing 
only a Ag-type two derivative term: 

ST = -cJ ds% [m E (<p) + N%«(<p) g*» v d, <pd v <p], (A4) 

where 

N*™ = N E + a E L E + 2I E . (A5) 

The conclusion is that, modulo 0(v 4 /c 4 ) terms and higher-order terms in the radius of the 
extended body E (associated to higher-derivatives), there is only one relativistic form factor 
for, possibly compact, extended bodies in tensor-scalar gravity: N E ew . By dimensional anal- 
ysis [iVg ew ] = [mass] [length] 2 and we therefore expects N E ew to be, roughly, some spherical 



inertia moment. Using the diagrammatic method of Ref. [12|, it is easy to compute the extra 



contribution to the A^-body Lagrangian entailed by the presence of A^: 

extended r ST G * ™A mg ^ A N ™ W Q 'b) n EA [ ^-EB ( 

A^E^B ' EA< EB C 

where hea = (%e — ^a)/^ea- Note that in a two-body system the summation will have two 
terms: one where E — 1, A — B — 2, and one where E — 2, A — B — 1. [We denote here 
for clarity the two bodies by the labels 1 and 2.] 

The above considerations have the advantage of being valid even when discussing strongly 
self-gravitating extended bodies. In the case of a weakly self-gravitating extended body, 
Nordtvedt |36| has obtained an effective action, after some explicit calculations of the equa- 



tions of motion, of the form 

1 



S E — — c I ds* E 



m E {y)+-I E A\ v )u tl u v R l 
o 



(A7) 



where Ie = J d 3 x p(x) x 2 is the physical spherical inertia moment, and where we recall 
that the tilde refers to physical, Jordan-frame quantities: g^ u = A 2 g*, R^ u = R pv [g], 
g flu u fJj u u = — 1. By expanding ( |A7| ) in terms of our generic expansion ( |A2| ) we find 

I E = , J E = — I E , K E = — r I E , L E = — Ie 
oo 6 

M E = -~ (A) - «o) Ie , Nf = ~ % + 2a 2 ) I E . (A8) 
6 o 

We conclude from Eq. ( |A~5| ) that the only observable form factor from infinity is 

NT= l -^I E . (A9) 

In other words, after the shift z' E = z E — | «o (Ie/^e) V* m if, the only new contribution 
to the A^-body Lagrangian is Eq. ( |A6| ) with N E ew given by ( |A9| ) . When comparing with 
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Nordtvedt's results, note that he does not perform the simplifying shift that we are advo- 
cating so that he ends up with a more complicated-looking iV-body Lagrangian containing 
a contribution like ( |A6| ) but with N B d = ~ (/3 + 2a 2 ,) I e , plus two other sums which are 
simply obtainable by applying the inverse shift z N°rdtvedt _ z us _|_ i q:o (-^/^.e) Vy>, with 
ip = fo — J2a m A ola/ta c 2 , to the unperturbed Lagrangian Lq(z A) va) = 12 a \ m A + 
J2a^b \ G 1f m A m B (1 + ot A «b)/^b- [Here, both v A = z A and r AB = \z A - z B \ get modi- 
fied by the shift.] Finally, we find that in a binary system the finite-size effects cause the 
appearance of the additional interaction energy terms 

^extended 771 1 ^^A^B , 9 a rnew 1 „ 2 a rnew \ /*in\ 
Ant = H 4 2 ("A^S + a B N A )■ ( A1( J) 

For compact bodies (radius R A comparable to G^m^/c 2 ) we expect that N™ w ~ 
/3A m A(G* fnA/c 2 ) 2 , so that the new interaction energy ( |A10| ) will be of the form 

^extended ^ _ ^ ^ ft 2 + ^ ^ ft 2 ) ^GR ? (An) 

where E^^ denotes the general relativistic third post-Keplerian (0(t> 6 /c 6 )) contribution to 
the interaction energy, and where C A and C B are numerical coefficients which are roughly 
of order unity. The 3PK energy ~ G * m A m B jr AB x (G*m A /c 2 r AB ) 3 depends on the 

distance r AB in the same way as S extended E int . 

The conclusion is that scalar-mediated finite-size effects can be neglected in a matched- 
filter analysis of the phase evolution of inspiralling binaries, because they only modify by 
a factor ~ 1 + C (3 a 2 (which tends to 1 as a 2 — > 0) terms already present in the general 
relativistic phase evolution. As discussed for similar fractional corrections in Section III, 
they can be neglected compared to the non-general-relativistic scalar dipolar contribution 
to the phase evolution. 
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FIGURES 




FIG. 1. Region of the (ao,/3o) theory plane allowed by solar-system tests, binary-pulsar exper- 
iments, and future gravity-wave detections, in the case where nuclear matter is described by the 
polytrope (^-(^l). In view of the reflection symmetry qo — ► — ao, we plot only the upper half 
plane. The region allowed by solar-system tests is below the thin line labeled "1PN". The PSR 
0655+64 data constrain the values of «o an d /?o to be between the two solid lines. The regions 
allowed by the PSR 1913+16 and PSR 1534+12 tests lie respectively to the right of the bold line 
and of the the dashed line. The horn-shaped region at the top-left of the dashed line is removed if 
the observable P^ bs is taken into account for PSR 1534+12. Each of these curves determines the 
level x 2 = 1 f° r the corresponding test. We have also plotted the level \ 2 = 2 for PSR 1913+16 
to underline that the precise value of x 2 is not very significant in the region where binary-pulsar 
experiments are more constraining than solar-system tests {[3q < — 3). The regions excluded by the 
gravity- wave observation limit ( |3.10| ), with a signal-to-noise ratio S/N = 10, lie on the hatched sides 
of the curves labeled "LIGO/VIRGO" . The case of a 1.4m0-neutron-star and a 10 mQ-black-hole 
binary system is labeled "NS-BH", whereas the case of a 1913+16-like binary-neutron-star system 
is labeled "NS-NS". The region simultaneously allowed by all the tests is shaded. 
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FIG. 2. Same plot as Fig. 1 in the case of a soft equation of state (Pandharipande) . The region 
possibly excluded by the LIGO/VIRGO detection of a (1.4m ) neutron star-(lOm ) black hole 
system lies above the dotted line. The bubble-like region at the left of Fig. 1 (binary-neutron-star 
system detected by LIGO/VIRGO) does not exist in the case of this soft equation of state. 




FIG. 3. Same plot as Fig. 2 for a medium equation of state (Wiringa et al ). Dotted lines 
indicate the regions excluded by future gravitational-wave observations, respectively inside the 
bubble for the NS-NS case, and above the straight line for the NS-BH case. 
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FIG. 4. Same plot as Fig. 3 for a stiff equation of state (Haensel et al). Note that the bubble 
excluded by the detection of a binary-neutron-star system by LIGO/VIRGO is much larger when 
the equation of state is stiff. 




FIG. 5. Same plot as Fig. 1, assuming the same polytropic equation of state, but a sig- 
nal-to-noise ratio S/N = 100 for the LIGO (hatched) curves. For clarity, the dashed line cor- 
responding to the PSR 1534+12 test has been suppressed. 
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